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ABSTRACT 


This  report  discusses  the  theoretical  impossibility  of 
deducing  the  surface  reflectivity  of  a rotating  illuminated  body 
from  its  lightcurve.  We  first  derive  the  lightcurve  for  an 
arbitrary  convex  body  as  an  expansion  in  spherical  harmonics. 
From  the  form  of  the  resulting  lightcurve  and  its  Fourier  series 
the  negative  result  follows  immediately. 
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I . INTRODUCTION 

The  analysis  of  the  lightcurve  of  a rotating  body  is  a 
topic  of  long  standing  in  astronomy.  The  natural  satellites 
of  the  other  planets  and  the  minor  planets  have  been  investigated 
in  this  fashion.  The  primary  aim  has  been  to  deduce  the 
surface  reflectivity  from  the  lightcurve  so  that  a deeper 
analysis  of  the  object  may  be  successful.  There  is  a result 

It 

due  to  Russell  concerning  the  impossibility  of  deducing  the 
surface  reflectivity  from  the  lightcurve.  For  the  natural 
satellites  this  constraint  is  of  little  importance  because 
the  most  sophisticated  model  would  be  an  airless  triaxial 
ellipsoid.  Data  acquired  by  spacecraft  have  shown  that  for 
some  moons  this  is  not  true,  but  then  we  have  the  flyby  data 
to  aid  in  our  analysis. 

In  the  case  of  artificial  satellites,  particularly 
those  covered  by  solar  cells,  high  order  symmetries  may 
exist  in  the  surface  reflectivity.  As  we  shall  see  below 
these  may  be  impossible  to  discover  from  only  the  lightcurve. 

We  expand  on  Russell's  result  and  fill  in  some  gaps  in  his 
analysis.  We  also  derive  the  Fourier  series  (in  time)  of 
the  lightcurve. 

*H.  N.  Russell,  As trophy  FT  J.  24,  1,  (1906). 
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II.  A SPHERICAL  BODY 

Let  r,  6,  4>  be  a spherical  coordinate  system  with  origin 
at  the  center  of  the  sphere:  r is  the  distance  from  the  origin, 

0 is  the  colatitude,  and  <J>  is  the  azimuth.  Let  B(0,4>)  be  the 

* 

intrinsic  brightness  of  an  element  of  area  dS(6,4>)  of  the 
surface.  Let  y(9,4>)  be  the  angle  between  the  outward  normal 
to  dS  and  the  line  of  sight  to  the  observer.  Finally,  let 
R(6  ,<(>)  be  the  distance  between  dS  and  the  observer.  The  amount 
of  light  received  by  the  observer  from  dS  is  its  brightness 
multiplied  by  its  apparent  area  divided  by  the  square  of  its 
distance;  viz.. 


dL  = BcosydS/R 


The  total  light  received  by  the  observer  is  the  sum  of  all  such 
contributions  for  which  cosy>o  (e.g.,  that  dS  is  visible). 


L = /BcosyH (cosy) dS/R  , 


where  H(u)  is  the  Heaviside  function 


H (u)  = 


1 u>o 
0 u<o' 


The  magnitude  is  related  to  the  brightness  by  m = -2.5£og(B) 
+ constant. 
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If  the  radius  of  the  sphere  is  a and  a<<min[R(0 , 4>)  ] then 

0 /<J> 

L = (1/R2) /BcosyH (cosy)dS,  (3) 

where  R is  now  the  distance  from  the  center  of  the  sphere  to  the 
observer  and  cosy  = r . R / (rR)  . The  vector  R = (R,0,4>) 
where  0 is  the  observer's  colatitude  and  $ is  his  azimuth 
(the  vertex  of  y is  now  at  the  origin) . Since 


and 


cosy  = cos0cos 0 + sin0sin0cos  (<}>-$)  , 
dS(0,$)  = a2sin0d0d()) , 


(4) 

(5) 


L = (a  /R  j /Qd0/oB  (0  , ip)  cosyH  (cosy)  sin0d<f> . (6) 

2 2 
Clearly  L = a (function  of  0,$>)/R  . 

To  pursue  the  analysis  further  we  now  assume  that  B(0,<(>)  is 
non-negative,  everywhere  bounded,  has  a finite  number  of 
extrema,  and  a finite  number  of  points  or  curves  of  discontinuity. 
These  conditions  (known  as  Dirichlet  conditions)  pose  no 
practical  limitations  on  the  form  of  B.  They  do  assure  us 
that  the  series  expansion  for  B,  in  a series  of  spherical 
harmonics,  converges.  For  reference  we  define  the  spherical 
harmonics  { Y£m C © » J in  terms  of  the  associated  Legendre 
functions  (P™(cos0)}  and  exponentials  (exp(im<(>)}  by 
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d 


\ 


we 


r -.1/2 

j.  v |2 1 + 1 r (t  — m + 1)  Lm  / _ ns  _/• 

= [-4?—  r'T-TT  m'  + T)  jP£  (cos0)  exp  (im<|))  , 


Pm  (U>  = ~ u2)m/2  + m(u2  - 1)£ 

1 2*TU  + 1)  " du£  + m 


(7a) 


(7b) 


= (-l)m(l  - u2)m/2dmP£(u)/dura, 


(7c) 


P"m(u)  = (~l)luT[l  - m + 1)  p“l(u)/rU  + m + 1)  , 


,mr 


,m, 


(7d) 


P™  (-u)  = (-1)1  + ^(u)  , 


(7e) 


??(1)  = 6 
£ mo 


■t 


1 if  m = 0 
0 otherwise  ' 


(7f ) 


Yl  - m(9  ,<^)  = 


Im 


(7g) 


where  r (u)  is  the  gamma  function  and  P^  (u)  is  the  Legendre 
polynomial  of  order  l. 

We  write  B (0,<f>)  as 


1=0  m=-l 


(8a) 


where 


Bfm  = /od0/o7rsin0Y*m(0,4>)B(0,4.)d<t., 


Im' 


(8b) 


and  the  asterisk  denotes  complex  conjugation.  We  substitute 
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r 


expression  (8a)  into  Eq.  (6)  and  deal  with  the  Heaviside  function 
by  expressing  in  terms  of  y and  A.  The  auxiliary  angle 

A is  the  position  angle  of  dS,  relative  to  the  origin,  as  seen  by 
the  observer  after  the  surface  is  projected  onto  the  plane  of  the 
sky.  From  the  properties  of  the  rotation  group  of  the  spherical 
harmonics  we  can  write 

l 

Y£I»(e"M  = nLtAinYen(v'M-  (9a) 

However,  from  Eq.  (4)  , if  9 = 0 and  <J>  = <t> , y = 0.  From  Eq.  (7f) 
it  now  follows  that 

Yfm(0'$)  = Afo-  (9b) 

We  now  use  Eqs.  (6,  8a,  9a,  and  9b)  to  write  the  f'th  order 
term  for  L as  [cf.  Eq.  (11)] 

Lf  “ <a2/R2)  l B tJ0  dA//  sinYCOSYlY£ra(0,*)Y£o(Y,A) 
m=-t 

* I,  W1'1111  •U14'-  (10a) 

n=— L 

From  Eq.  (7a)  we  see  that  the  result  of  the  A intergral  is 

2tt6  ^ so 
no 
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(10b) 


I<£  * (2Tra2/R2)  £ B^mY^m(0,4>)/^/zsinYCosYY£o(Y,X)<3Y, 
m=-£ 


r"/2. 


where 


LZo 


[2 1 + iii/2„o,  x r 21  + 1]  i//2„ , 

= [-TF-J  Vcosy)  =[“4i J p£  (cosY) 


If  we  set  I£  equal  'to  the  integral  in  Eq.  (10b)  then 

l 

— 2 T /n2x 


L,  = (2,a2VR2)  £ BtmYta(0,t), 
m=— -t 


(10c) 


(10d) 


and. 


OO  00  £ 

L = I L.  = ( 2ira2/R2)  I I£  I B/mYto ( 0 » * > 
£=o  * l=o  m=-l  *** 


(11) 


Equation  (11)  is  our  final  result  giving  L as  a function  of  the 
(B£m>  and  the  observer's  location  R. 

The  last  step  is  the  evaluation  of  the  integrals  {l£}.  This 
is  done  by  using  the  Rodrigues  formula  for  the  Legendre  poly- 


1 


4tt 


SL/2 


21+1 


1l  ~ /^uP^fuidu  = 


1/2  if  l = 0 
1/3  if  i = 1 
1/8  if  l = 2 


(12) 


0 if  Z is  odd  and  >3 

1 

2 + 1 


(~1-)- (-£-+^yiy'li  if  1 is even  and ^ ' 


where  the  double  factorial  is  defined  by 


(2k)  ! ! = 2*T  (k  + 1) 


(2k  + 1)!!  = F(2k  + 2) / [ 2 T (k  + 1)] 


k = 0,  1,  2, 


(13) 


There  are  two  conclusions  we  can  immediately  draw  from 
Eq.  (12);  first  for  non-zero  values  of  I, 


1 IU  + 2)  1 „2  V*-2' 


(14) 


so,  since  the  series  for  B [Eq.  (8a)]  converges,  the  series  for 
L [Eq.  (11)]  converges  uniformly  and  absolutely.  Hence,  L is  a 
continuous  function  of  R and  V^L  is  continuous  even  if  B is  dis- 
continuous. Secondly,  since  all  of  the  higher  order  odd  harmonics 
are  absent,  there  exists  a denumerable  set  of  different  reflec- 
tivities all  of  which  yield  the  identical  light  curve. 
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III.  EXTENSION  TO  ANY  CONVEX  SURFACE 


The  simplest  two-dimensional  surface  to  discuss  is  one 
whose  equation  is  given  in  the  form  z = f(x,y).  However,  we 
will  use  the  parametric  representation  r = [x(0,<}>)»  y ( 0 , 4> ) , 
z (0  ,<J>)  ] . We  assume  (i)  r(0,4>)  is  a single  valued  function, 
(ii)  dr/dQ  and  3r/d<f>  both  exist  and  are  continuous,  and  (iii) 
dr  dr 

-rs-  x is  not  null  (i.e.,  the  vector  cross  product), 

dt)  09 

The  element  of  arc-length,  ds , on  the  surface  is  given  by 


dr 


dr 


or. 


ds  = dr  ' dr,  dr  = d0  + dtp , 


ds2  = Edd 2 + 2Fd0d <p  + Gd<p2 , 


where 


3£ 

E = d0 


0£ 

F = 

30 


dr 


30  ' 

ar 
dj 

ar 


(=  0 here) , 


G 34>  3<t>  * 


(15a) 


(15b) 


(16a) 


(16b) 


(16c) 


The  element  of  area,  dS,  of  the  surface  is  proportional  to 




I 3r  9ri 

aT  x 3^  ’ since»  for  any  two  vectors  u and  v, 

|u  x v|2  = |u|2|v|2  - (u  * v)2,  (17) 

dS  can  be  written  as 

dS  = (EG-F2)1/2d9d(J>.  (18) 

Finally,  define  a quantity  C (which  is  essentially  a curvature) 
by 

C = r2sin0/ (EG  - F2)1/2.  (19) 

Then  the  expression  for  L is 

L = J277/^  (B/C)  cosyH  (cosy)  sin0d<f>d0/R2 , (20) 

o o 

2 2 

in  analogy  with  Eq.  (lb).  If  max(|r|  )<<  min  (|r|  ) then 

9 / 4*  -j  2 9/4* 

we  may  again  bring  the  factor  of  a /R  outside  of  the  integral, 
redefining  R and  y as  in  8 II,  and  proceed  exactly  as  before 
except  the  {B^}  are  now  interpreted  as  the  expansion  coefficients 
of  B/C.  Only  if  C vanishes  will  our  previous  results  be 
invalid.  The  restriction  to  convex  surfaces  precludes  this 
possibility. 

We  might  inquire  as  to  the  lightcurve  of  an  arbitrary 
convex  body  with  a constant  surface  brightness.  From  the 
discussion  above  the  answer  is  again  given  by  Eq.  (11) 

(after  reinterpretating  the  { B^}  ) except  there  can  be  no  t * 1 
term.  This  follows  from  the  parity  of  the  spherical  harmonics, 

Eq.  (7g) . 
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IV.  THE  FOURIER  SERIES  FOR  THE  LIGHTCURVE 


As  the  orientation  of  the  r,0,<|>  coordinate  system  is  at  our 
disposal  we  now  choose  the  polar  axis  to  be  the  axis  of  uniform 
rotation.  If  the  body  is  precessing  then  the  polar  axis  is  the 
instantaneous  axis  of  rotation.  The  observer's  colatitude  is 
a constant  equal  to  the  inclination  of  the  axis  of  rotation  to 
the  line  of  sight.  The  observer's  azimuth  is  simply 
4>  = a)  ( t - tQ)  where  10  is  the  angular  speed  of  rotation  and  tQ  is 
the  instant  of  a prime  meridian  crossing.  Using  this  in  Eq.  (11) 
and  rearranging  terms  (which  is  permissible  because  the  series 
converges  absolutely)  we  obtain 


L = (2ira2/R2)  (C  (0)/2  + l C (0)cos[nu>(t  - t )] 


n=l 


n 


+ l S (0)  sin  Inu)  (t  - tQ)]}, 
n=l 


(21) 


with 


From  Eqs. 


co(e)  ■ \l  'AA1™81' 

K“0 

cn  ' j V'Wk'0030’  + VA”*00*011-  n>° 
k=n 

sn  “ J iIk[BknIk(cos9>  ' Bk-npk"(cos0)1'  n>°- 
k=n 

(?g,  8b),  Bk_n  = (-l)Xn. 


(22a) 


(22b) 


(22c) 
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V.  THE  INVERSE  PROBLEM 


As  we  mentioned  above  there  exist  an  infinite  number  of 
different  B's  all  of  which  yield  the  same  L.  Clearly  then, 
given  L it  is  impossible  to  determine  B uniquely.  Any 
additional  physical  constraints  we  might  want  to  impose  on  B 
(i.e.,  its  non-negativity,  continuity,  various  symmetry 
properties,  etc.)  can  serve  to  limit  our  ignorance  but  not 
overcome  it.  The  problem  is  simply  indeterminate.  As  this 
is  the  case  when  a complete  set  of  functions  is  used  the  situation 
will  be  even  worse  if  some  other,  incomplete  basis  is  used. 
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